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Abstract 

Particle production at the end of a first-order electroweak phase transition may 
be rather generic in theories beyond the standard model. Dark matter may then be 
abundantly produced by this mechanism if it has a sizable coupling to the Higgs field. 
For an electroweak phase transition occuring at a temperature Tew ~ 50 — 100 GeV, 
non-thermally generated dark matter with mass Mx > TeV will survive thermalization 
after the phase transition, and could then potentially account for the observed dark 
matter relic density in scenarios where a thermal dark matter component is either 
too small or absent. Dark matter in these scenarios could then either be multi-TeV 
WIMPs whose relic abundace is mostly generated at the electroweak phase transition, 
or "Baby-Zillas" with mass Mqut ^ Mx ^ ^ew that never reach thermal equilibrium 
in the early universe. 



1 Introduction 



The most popular paradigm for the origin of dark matter (DM) in the Universe is the 
thermal freeze-out. In that scenario, the dark matter particle with mass Mx annihilates 
into matter with a cross section (a v ) thcrmal ~ 3 x 10~ 26 cm 3 /s. This ensures dark matter 
is in thermal equilibrium with the rest of the plasma in the early universe while T > M 
but decouples when T ~ Mx/20, leaving the relic abundance in agreement with the value 
Qx = 0.228±0.027 measured by WMAP [I]. Incidentally, (a w) thcrmal is a generic cross section 
for a weak scale mass particle interacting with order one couplings, this fact being referred to 
as the WIMP miracle. In spite of these attractive features, non-thermal mechanisms of dark 
matter production have also received considerable attention. Examples include right-handed 
neutrinos produced by oscillations [2], axions produced by vacuum misalignment [3], winos 
produced from moduli decays [I], and super-massive dark matter (WIMP-zillas) produced 
during reheating after inflation [5]. These studies allow one to recognize a wider range of 
possible collider and astrophysical signals of dark matter than what would result from the 
thermal WIMP scenario. 

In this paper we study the possibility of non-thermal dark matter production during a 
first-order electroweak (EW) phase transition. Bubble collisions at the end of the EW phase 
transition may give rise to abundant non-thermal particle production when a sizable amount 
of the energy budget of the transition is stored in the bubble walls, possibly leading to new 
and appealing scenarios. Many models of dark matter contain a direct coupling between the 
Higgs and the dark matter candidate fields (MSSM and its extensions, Little Higgs theories 
with T-parity and DM extensions of the standard model (SM) via the Higgs portal, to 
name a few). It is thus reasonable to expect that dark matter may be abundantly produced 
non-thermally at the end of a first-order EW phase transition. Note that, much like in the 
thermal WIMP case, dark matter would then be a particle with Mx ~ 10 GeV — 10 TeV with 
significant coupling to the SM, thus being within reach of colliders and DM direct detection 
experiments. 

There is however one generic problem with this scenario. Since the temperature of the 
Universe right after the EW phase transition is Tew ~ 50 — 100 GeV (for strong transitions 
Tew may be somewhat lower than 100 GeV), thermalization will typically lead to a wash- 
out of the non-thermal abundance, thus rendering the particle production at the EW phase 
transition irrelevant for the subsequent evolution of the Universe. The wash-out process 
can nevertheless be avoided in a number of ways, each resulting in a scenario where non- 
thermal dark matter production is (fully or partially) responsible for the observed dark 
matter relic density. One possibility, recently outlined in [5J, is to allow for a few e-foldings 
of inflation prior to the beginning of the transition (which can happen for very strong EW 
phase transitions), diluting the plasma and leaving the Universe partially empty. If the 
reheating temperature after the phase transition is low, Trh *C 100 GeV, it may be possible 
for a dark matter candidate with weak couplings to the Higgs field and mass Mx ~ 100 
GeV to remain out of thermal equilibrium after the EW phase transition. In this paper we 
investigate other scenarios allowing for a survival of the non-thermal abundance. 

One possibility corresponds to the case of relatively heavy (multi-TeV) dark matter: for 
Mx > 1 TeV, dark matter will be very non-relativistic when the EW phase transition takes 
place, and the decoupling/freeze-out temperature T FO will satisfy T FO ~ Mx/20 > T E w- 
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Then, heavy dark matter produced non-thermally through bubble collisions may remain out 
of thermal equilibrium after the EW phase transition (or at least wash-out will be partially 
avoided). Another possibility is that bubble collisions produce super- heavy dark matter, 
Mx ~ 10 6 -10 8 GeV, a scenario we call "baby-zillas". We argue this may be possible for a 
very strong EW phase transition and dark matter having a large coupling the Higgs. In 
order for baby-zillas with Mx 3> ^ew to be a viable dark matter candidate, they must 
have never reached thermal equilibrium in the early universe after inflation, since otherwise 
they would have over-closed the universe. This sets a relatively low upper bound on the 
reheating temperature after inflation in that scenario. Finally, asymmetric dark matter 
production might allow to avoid complete wash-out of the non-thermal abundance through 
thermalization after the EW phase transition. 

The paper is organized as follows: in Section 2 we review the formalism that describes 
particle production at the end of the EW phase transition for the case of very elastic bubble 
collisions [TIE] and extend it to the case of very inelastic ones, highlighting the differences 
between both scenarios [9]. Then, in Section 3 we explicitly compute the particle production 
efficiency of scalar, fermion, and vector boson particles coupled to the Higgs (either directly 
or via an effective Higgs portal). In Sections 4 and 5 we focus on dark matter production 
at the end of the EW phase transition. First we discuss in Section 4 the conditions for 
non-thermally produced dark matter to avoid subsequent wash-out and constitute the bulk 
of the present dark matter density, selecting heavy (multi-TeV) vector boson dark matter as 
a viable example. We go on to analyze in detail non-thermal dark matter production in that 
scenario and the subsequent evolution of the non-thermally generated abundance after the 
EW phase transition, including finally a discussion on the current XENON100 bounds and 
direct detection prospects. Then, in Section 5 we study the non-thermal production of very 
heavy (Mx 3> vew) vector boson dark matter, and outline the conditions under which these 
baby-zillas constitute a viable dark matter candidate. In the case of asymmetric non-thermal 
dark matter production, we find it difficult to avoid subsequent wash-out, and the discussion 
is left for an appendix. We summarize our results and conclude in Section 6. 

2 Particle Production at the EW Phase Transition 

2.1 Bubble Collisions in the EW Phase Transition 

If the early Universe was hotter than Tew ~ 100 GeV it must have undergone an EW phase 
transition at some point in its history. Within the SM, the EW phase transition is a smooth 
cross-over, however it is conceivable that new degrees of freedom beyond the SM modify the 
Higgs potential so as to make the transition first order. This is what we assume throughout 
this paper, without specifying the full theory that makes the first order transition possible. In 
that case, the EW phase transition proceeded through nucleation and expansion of bubbles 
of true Higgs vacuum, which eventually collided among each other completing the transition. 
As this was happening during the radiation dominated era, the bubble expansion process 
would then have taken place in a thermal environment (except for the case when a period 
of inflation would have preceded the phase transition). 

For a first order phase transition occuring in a thermal environment, the study of the 
bubble expansion process reveals that the thermal plasma exerts some amount of friction 
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on the expanding bubble wall, and this friction tends to balance the pressure difference 
on the bubble wall driving the expansion. In the usual picture, nucleated bubbles reach 
a stationary state after a very short period of acceleration, with a constant wall velocity 
depending specifically on the interactions of the bubble wall with the degrees of freedom in 
the plasma [10|[TTj and on the resulting fluid dynamics [T2HH] (see [15] for a review). In this 
case, the amount of energy stored in the bubble walls at the time of the bubble collisions 
is negligible compared to the available energy of the transition, since most of this available 
energy gets converted into plasma bulk motion and thermal energy [16] . 

However, this picture was recently challenged in [17], where it was shown that the friction 
exerted by the plasma may saturate to a finite value for ultrarelativistic bubble walls. As 
a consequence the stationary state assumption will no longer be true when the pressure 
difference on the bubble wall is larger than the friction saturation value, which may happen 
for strongly first order phase transitions. In this scenario, if there are no hydrodynamic 
obstacles that prohibit the bubble walls to become highly relativistic in the first place (see 
however [IE]), bubbles will expand in an accelerated way ('the so-called runaway bubbles), 
with almost all the energy of the transition being used to accelerate the bubble wall£] [T5] . 
By the end of the phase transition (when bubbles start colliding), these runaway bubbles 
may reach very large values of 7^: 

ry,„ < 7 max ~ ^- MbL rsj 10 15 (2 1) 

with vt the value of the Higgs VEV in the broken phase and ~ (10~ 3 — 10 -2 ) if -1 being 
the duration of the phase transition |19j. The estimate (12. ip follows from balancing the 
surface energy on the bubble wall ( 12. 5 j) and the energy available inside the bubble. 

Once bubbles start colliding, the energy stored on the bubble walls will be liberated into 
the plasma. As argued above, for "runaway" bubbles this will correspond to a very large 
portion of the energy budget of the phase transition, and therefore this process can be very 
important. Under certain circumstances, this may also hold true for highly relativistic bubble 
walls (tu, ^> 1) that reach a stationary state long before bubble collisions start (meaning 
that 7^ <C 7™ ax ), in which case the amount of energy stored in the bubble walls will be very 
small compared to the available energy of the transition, but still important when released 
into the plasma at the end of the transition. 

The process of bubble collisions in cosmological first order phase transitions is by itself 
a very complicated one. Consider a configuration of two planar bubble wall^l initially far 
away from each other, that approach and collide [7J[8j[20]. Depending of the shape of the 
potential for the scalar field <p driving the transition (in our case, the Higgs field h), the 
bubble collision will be approximately elastic or partially inelastic [TIEj (see also [9]). In the 
first case, the walls reflect off one another after the collision, which reestablishes a region of 



lr TIiis situation may also arise if, under very specific circumstances, a few e- foldings of inflation are 
achieved prior to the beginning of the EW phase transition (see [5] for a natural realization of this scenario), 
diluting the plasma and leaving the Universe mostly empty. In this case the expansion of the bubbles 
effectively takes place in vacuum, and the nucleated bubbles expand in an accelerated way due to the 
absence of friction. 

2 At the time of the collision, the bubbles are so large compared to the relevant microscopical scales, that 
their walls may be considered planar as a good approximation. 
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symmetric phase between the bubble walls. For a perfectly elastic collision the field profile 
of the colliding walls in the limit of infinitely thin bubble walls (taken as step-functions) can 
be written as [8] 

{0 if v w t < z < —v w t t < 0, 
if - v w t < z < v w t t>0, (2.2) 
Vt Otherwise, 

where v w is the bubble wall velocity, the bubble walls move in the ^-direction and the collision 
is assumed to happen at t — 0. Since we are ultimately interested in scenarios where j w 3> 1, 
we will take the ultrarelativistic limit v w — > 1 in the rest of the section. The field profile 
(12. 2p neglects the thickness of the bubble walls l w (generically, l w ~ (10 — 30)/T EW , with 
Tew ~ 50 — 100 GeV). To capture the wall thickness effects one can consider another ansatz 
for the profile of the colliding bubble walls: 



h(z,t) = h lv 



Tanh ( 7, 



Vt 

y 



Vt 

T 



2 + Tanh 7 



t + Id 



z - \t\ 



Tanh ( 7 



Tanh ( 7 



t - \z\ 



z + \t\ 



(2.3) 



A perfectly elastic collision is however an idealized situation, as one expects a certain 
degree of inelasticity in a realistic collision. Moreover, even for a very elastic collision the 
bubble walls will eventually be drawn back together by vacuum pressure and collide again. A 
quantitative picture of the collision of two planar bubble walls can be obtained by studying 
the evolution equation for the scalar field configuration h(z, t) subject to the potential V(h): 



{%-%)h{z,t) 



dV(h) 



(2.4) 



with the initial condition corresponding to two planar bubble walls far away from each other 
and moving in opposite directions (given approximately by hi w in the limit t — > —00). In the 
ultrarelativistic limit the ansatz (I2.3P will also be an approximate solution of (12.41) before 
the bubble collision^! (for t < 0). In this limit, the kinetic energy per unit area contained in 
the field configuration h(z, t) prior to the collision is given by 



E w 



2 2 lw 

3 Vt ~- 



(2.5) 



At the moment of the collision, the field configuration makes an "excursion" to field values 
larger than vt in a small region around the collision point [20] (resulting in dVj dh 7^ in this 
region). The subsequent evolution of h(z,t) strongly depends on the shape of the potential 
V(h). The field close to the collision region oscillates back after the initial "kick" in field 
space, and for a potential with nearly degenerate minima this oscillation is able to drive the 
field over the potential barrier and into the basin of attraction of the symmetric minimum 



3 Each bubble wall in (|2.3p interpolates between the symmetric and broken minima of V(h), and so 
dV(h)/dh = outside the bubble wall. Then, for very thin walls the equation of motion approximately 
simplifies (before the collision) 
f(z — t) is a solution. 



to (<9 t 2 — d^) h(z,t) = 0, for which any function of the form f(z + t) or 
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Figure 1: LEFT: Potential with nearly degenerate minima. RIGHT: Potential with very non- 
degenerate minima. Each one shows the behaviour of the field immediately after the collision in 
the region close to the collision point, as described in the text: 1) "Kick" to field values larger 
than v(T). 2) Large field oscillation, successful (LEFT) or unsuccessful (RIGHT) in driving the 
field over the potential barrier. 3) Oscillations around the symmetric (LEFT) or broken (RIGHT) 
minimum. 



(Figured]- Left), where it will perform small-amplitude oscillations. In this case the collision 
is approximately elastic as described above, with the bubble walls being effectively reflected 
as a region of symmetric phase is re-established between them. The walls move then away 
from each other until vacuum pressure makes them approach and collide again, repeating 
the process several times. In each collision some fraction of the energy stored in the walls 
is radiated into scalar waves and quanta of the fields coupled to h, until all of the energy 
in the walls is radiated away. In contrast to this scenario, for a potential V(h) with very 
non-degenerate minima (Figured]- Right), the field oscillation after the "kick" in the region 
close to the collision point does not effectively drive the field over the potential barrier. As 
a consequence, the field stays in the basin of attraction of the broken minimum vt and 
performs relatively large- amplitude oscillations around it, giving rise to a large amount of 
energy radiated into scalar waves (as opposed to the previous scenario). In this case the 
collision is very inelastic. 

Following [9], we compute the numerical solution for the field profile h(z, t) corresponding 
to the collision of two bubble walls, obtained from solving ( 12. 4ft with a toy potential V(h) of 
the form 

V{h) = a 2 h 2 - b 2 h 3 + A/i 4 (2.6) 

both in the case of nearly degenerate minima (Figure 1 - Left) and very non-degenerate 
minima (Figure 1 - Right). The results are shown in Figure [2] (similar plots appeared 
earlier in [21]). Figure 2 - Left (corresponding to the potential of Figure 1 - Left) shows an 
approximately elastic bubble collision, while Figure 2 - Right (corresponding to the potential 
of Figure 1 - Right) shows a very inelastic one. 

Guided by the numerical solution for h(z, t) in the case of a very inelastic collision, we 
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h(z,t) 




h(z,t) 




Figure 2: Snapshots of the field profile h(z,t) during a bubble collision (t increasing downwards). 
LEFT: Bubble collision for the potential with nearly degenerate minima (Figure 1 - Left). RIGHT: 
Bubble collision for the potential with very non-degenerate minima (Figure 1 - Right). In both 
cases, j w = 10 2 , l w = 15/Tew and Tew = 100 GeV. 



can obtain an analytic solution h(z,t) = h^i for the case of a "totally inelastic collision" 
(as opposed to the "perfectly elastic collision" described earlier), in which all the energy is 
radiated in the form of scalar waves after the bubble collision. For t < (before the collision) 
we have 



h T i(z,t < 0) = v T + 



Vt 



Tanh ( 7 



z + 1 



Tanh ( 7 



t 



l I: 



(2.7) 



which matches hi w (z, t < 0). In order to obtain fiTi(z, t) for t > 0, we note that the field will 
not leave the basin of attraction of the broken minimum vt after the collision. We can then 
approximate the potential V(h) the field will feel for t > as 



V(h) 



mi 



(h - v T )' 



mi 



(2. 



2 " ' 2 

This allows to solve the equation of motion (12.41) explicitly for 8tiTi(z,t) = h^\{z,t) — vt- 



(d 2 t - dl) Shn{z,t) = -m 2 h Sh Tl (z,t) 



6h T i(z,0) = h lw (z,0)-v T = 



(2.9) 



d t 6h T i(z,0) 



Cosh 2^ 
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where the boundary conditions follow from imposing continuity of 8h^i(z, t) and d t 5h,Ti(z, t) 
at t — 0. From (12.91) . we finally obtain 



hTi(z, t > 0) = vt 



In 



dp z 



Pz 



Cos (p z z) 



VpI + m l Sinh ( ^ 



Sin 



(2.10) 



Notice that in the limit m h — > 0, (I2.9p becomes (<9 2 — d 2 z ) <5/zti(^, t) = and (12. 7p is also a 
solution for t > 0, case in which the two bubble walls would pass through each other without 
actually colliding. 

The analysis for the dynamics of bubble collisions presented here may be extended to 
phase transitions involving multiple fields (see for example [2D]), although in this case the 
analysis of the field evolution after the bubble collision becomes much more complicated 
(since the scalar potential is multidimensional and the field "excursion" at the moment of 
the bubble collision will involve several fields), and we will not attempt it here. 



2.2 Particle Production Through Bubble Collisions 

The bubble collision processes analyzed in the previous section allow to liberate into the 
plasma the energy contained in the bubble walls. This can happen either via direct particle 
production in the collisions or via radiation of classical scalar waves which will subsequently 
decay into particles. For bubble collisions taking place in a thermal environment, the number 
densities n a for the different particle species created during the collisions should very quickly 
approach the ones in thermal equilibrium after the phase transition, thus rendering the 
particle production process irrelevant for the subsequent evolution of the Universe. However, 
as it has been briefly discussed in the introduction, under certain conditions fast thermal- 
ization of certain species after the phase transition may be avoided, which can make the 
particle production process very important in that case. 

In order to study the particle production through bubble collisions, we will treat the scalar 
field configuration h(z, t) as a classical external field and the states coupled to it as quantum 
fields in the presence of this source. In doing so, we will neglect the back-reaction of particle 
production on the evolution of the bubble walls themselves throughout the collision, which 
should be a good approximation when the energy of the produced particles (for each species) 
is much less than the energy contained in the field configuration h(z,t). The probability of 
particle production is given by [8] 

V = 2lm(T[h}) (P<1) (2.11) 

where T [h] is the effective action. T [h] is the generating functional of 1PI Green functions, 
and to the quadratic order in h 

T[h] = ^J d 4 Xl d 4 x 2 h( Xl ) h(x 2 ) T (2) (xi, x 2 ) (2.12) 

with r 1 - 2 ) (xi,x 2 ) = r 1 - 2 ) (xi — x 2 ) being the 2-point 1PI Green function. In terms of its 
Fourier transform T^ 2 ^ (p 2 ), and using (12. lip and (12. 1 2 j) we get 
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V 



j T^Im (f (2) (p 2 )) J d 4 Xl d 4 x 2 h{ Xl ) h(x 2 ) e ip ^~ X2) (2.13) 



The last integral in ( 12. 13ft is just h(p) , with h(p) being the Fourier transform of the 
Higgs field configuration h(x) 



Hp) 



d 4 xh(x) e ipQ 



(2.14) 



For a background field configuration h(z,t), its Fourier transform is given by h(p) = 
(27r) 2 S(p x ) S(p y ) h(p z ,cu). Then, using (I2.13p . we obtain the mean number of particles pro- 
duced per unit area [8]: 



dp z dui 



h(p z: u)\ 2 lm (f( 2 ) (u 2 -p 2 z ) 



(2.15) 



The physical interpretation of (I2.15P is rather simple [8]: the scalar field configuration 
h(z, t), corresponding to the two bubble walls that approach and collide, can be decomposed 



into modes of definite four-momentum p 



PI via the Fourier transform. Modes 



with p 2 > represent propagating field quanta with mass squared m 2 = p 2 . Then, (12.151) 
integrates over the amount of field quanta of mass p 2 contained in the field configuration 
multiplied by the probability of those quanta to decay. 

The Fourier transform of the background field configuration h(z, t) can be performed 
explicitly both for the case of a perfectly elastic collision and of a totally inelastic one 
analyzed in the previous section. For a perfectly elastic collision, in the limit of infinitely 
thin walls (h(z,t) = h^), we obtain 



h(p z ,u) = h oc (p z ,u) = 



Pz 



(2.16) 



However, since the highest values of p z and u available in the field configuration are 
naively expected to be of order 7 tu // iu (modes with p z ,u 3> ^y w /l w will be exponentially 
damped), the integration in (I2.15P should in this case be cut-off for p z > j w /l w and uj > 
Jw/lw From (I2.15P and (I2.16P we then obtain 
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(2.17) 



{u 2 -p 2 z ) 2 

Alternatively, when the thickness of the bubble walls is accounted for (h(z, t) = hi w ), the 
Fourier transform of (12. 3 p gives 



h(p z ,u) = h lw (p z ,u) 



re L, u 



2 ^ Sinh 



n l w lu 
2"/ w 



ui 2 — p 2 z 



(2.18) 



which automatically incorporates the exponential damping for u,p z 3> J w /lw The mean 
number of particles per unit area now reads 
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(2.19) 



For the opposite case of a totally inelastic collision (h(z,t) = /in), the Fourier transform 
is given by 



h(p z ,uj) = h T i(p z ,uj) = 



2 7™ Sinh 



2v T 



27™ 



uj 2 -p 2 z 



uj 2 -p 2 z 



m 



(2.20) 



The relative "— " sign between the two contributions in (I2.20p can be easily understood 
noticing that in the limit m h — > the Fourier transform of h T i(z,t) should give h(p z ,uj) ~ 
5(uj ±p z )- From (I2.20p . the mean number of particles produced per unit area in the case of 
a totally inelastic collision is given by 
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T lw Pz 
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(2.21) 



The expressions (12.171) . (12.191) and (12.211) can be rewritten in a more compact form by 
making the change of variables X = ^ ~ V%i ^ = ^ + v\- After performing the integral in 
\&, the mean number of particles produced per unit area finally reads 



jV 



A 2tt 2 



d X f(x)lm (f(2) ( x ) 



(2.22) 



The function f{ X ) encodes the details of the bubble collision process and quantifies the 
efficiency of particle production. For a perfectly elastic collision, in the limit of infinitely 
thin bubble walls, we have 



16 v\ Log 



fix) = foo(x) 



2fe) 2 -X+2^V / (e) 2 -* 
X 
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T J -x 



(2.23) 



For a perfectly elastic collision, and for bubble walls with finite thickness, we have 

* + X 1 



2tt 2 1 2 v 2 1 



Iw A J \ 
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2 v^7» 



(2.24) 



Finally, for a totally inelastic collision, we have 



fix) = hiix) = 



2 Tw X 2 ix-m 2 h ) 
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* 2 - X 2 ( Sinh 
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2v^7» 



(2.25) 
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In Figure [3] we compare the efficiency f(x) for the various cases (12.231) . (12.241) and (I2.25p . 
Notice that /ti(x) diverges as x ~ ► m \- This divergence is artificial, due to considering 
h(z,t) over infinite time and space, and should be cut-off since our solution is not valid 
over distances larger than the bubble radius Rb- Implementing this cut-off can be well 
approximated by replacing in ( 12.241) 



(x-ml) 2 -> (x-m 2 h ) 2 + {mill) h 2 w . 



(2.26) 




Figure 3: Particle production efficiency /(% = oj 2 — p 2 ) for j w = 10 2 (LEFT) and j w = 10 3 
(RIGHT), l w = 15/Tew an d Tew = 100 GeV, in the case of a perfectly elastic collision with 
infinitely thin bubble walls (|2.23p (solid red) and with a finite bubble wall thickness (|2.24p (dashed- 
black), and in the case of a totally inelastic collision (|2.25p (solid blue) with = 125 GeV. The 
X-axis is displayed in units of (100 GeV) 2 . 



Defining Xmm as the minimum value of x f° r which particle production is possible (cor- 
responding to the squared sum of the masses M a of the particles being produced), we im- 
mediately see from Figure E] that for a totally inelastic collision, production of light particles 
(Xmin < m D m ay be very efficient, while production of heavy particles (Xmin 3> m D wm 
be extremely suppressed. For a perfectly elastic collision, however, the production of heavy 
particles may be relatively efficient (we will comment further on this point at the end of 
section [3]). For the study of the efficiency of particle production in varios different scenar- 
ios in the next sections, we will use ( I2.23P for the case of an elastic collision, while for the 
case of a very inelastic one it is possible to show that ( I2.25P (together with ( I2.26P ) can be 
approximated as 



/ti(x) ^^v 2 T m\ 



Log 
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fe) 2 +x+ 2 evfe)'+x" 




X 




X 2 


(x-m 2 h f + m£4 

Iw 





(2.27) 



Let us now turn to the evaluation of the imaginary part of the 2-point 1PI Green func- 
tion's Fourier transform (% = cu 2 — p 2 ). Through the optical theorem, we can write: 
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Im (f( 2 ) (X))=\Y1 J dK a \M{h 



-> a) I 6 [x - Xr, 



(2.28) 



|2 . 



where .M(/i — >• a)\ is the spin-averaged squared amplitude for the decay of h into a set of 
particles a with masses M a , Xmin = (X) -^a) i s the minimum value of x for which this decay 
is possible and dH a is the relativistically invariant n-body phase space element 

d 3 h 1 



dli r 



n 



(2 7T) 3 2Ei 



(2 7r) 4 5 4 (p-^A; i ) 



(2.29) 



Then, the number of particles of a certain type a produced per unit area during the 
bubble collision will simply read from (12. 22ft and ( 12. 28 ft 



A 



4vr s 



Xmin 



dxfix) / rfn « \M(h ->■ a) 



(2.30) 



The amount of energy produced per unit area in the form of particles a is obtained by 
weighting (I2.30p by the energy of each decaying Fourier mode. This yields 



£_ 

A 



1 



4vr s 



dxfMVx dU a \M{h^a) 



(2.31) 



From (I2.30p and (I2.3ip . the non-thermally produced energy density p a (assuming that 
the produced particles quickly diffuse into the bubble interior) reads 

(2.32) 

v a a a v i± a tn B 

with A ~ 4 7r R 2 B being the total collision area and V the volume of the two colliding bubbles. 
From (12.321) . and bearing in mind that Rb — the non-thermally generated comoving 
energy density is 
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s(T E w) v 7 ^* m pi t ew H A 
with s(T EW ) the entropy density after the EW phase transition. 



(2.33) 



3 Particle Production via the Higgs Portal 

The efficiency of particle production may strongly depend on the nature of the particles 
being produced. In this section we will analyze the particle production efficiency for scalars 
S, fermions / and vector bosons coupled to the Higgs field. Apart from estimating the 
production of SM fermions and gauge bosons through this process, we will consider a simple 
Higgs-portal extension of the SM in order to study the production of other possible scalar, 
fermion or vector boson particles. Furthermore, we will restrict ourselves to Z 2 symmetric 
Higgs-portal scenarios, since we will ultimately be interested in dark matter analyses. We also 
comment on how to interpret the results in the case when the calculated particle production 
exceeds the energy available in the bubble wall. 
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3.1 Scalars 



For the complex scalar S interacting with the SM via the Higgs portal, the relevant part of 
the lagrangian is given by 



- AC S = m 2 \S\ 2 + X S \H\ 2 \S\ 2 with H 





~7T 



(3.34) 



2 

In this case, \M(h SS) \ = A>|, and one immediately obtains 
f (2) (X) 



Im 



J s 



\ 2 2 

\ s v T 



dU, 



1-4- 



6( X -4M 2 ) 



(3.35) 



with M 2 = m 2 + (A s /2) u| being the scalar squared mass. Then, using (12T221 . f[2T23|) . fl2~27j) . 
f)2.33p and (I3.35P we can compute the S-scalar comoving energy density generated through 
the bubble collisions (normalized to the observed dark matter comoving energy density) as 
a function of M s and X s . The results are shown in Figure HI 
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Figure 4: Generated comoving energy density T in scalars S (normalized to the observed dark 
matter comoving energy density) as a function of the scalar mass M s in the perfectly elastic collision 
limit (LEFT) and totally inelastic collision limit (RIGHT) for y w = 10 8 , l w = 15/Tew an d Tew = 
100 GeV. The solid-black line corresponds to the observed dark matter comoving energy density, 
and the dashed-black line (not seen in LEFT) corresponds to the maximum possibly generated 
comoving energy density (£ = E w ). 



From Figure H] it can be clearly seen that scalar particle production is quite suppressed for 
elastic collisions. For very inelastic collisions, heavy-scalar particle production is extremely 
suppressed, while production of light scalars turns out to be very efficient in this case. In 
fact, Figure H] shows that for large values of A s (A s < 1) the naively calculated energy of the 
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produced particles 8 exceeds the amount of energy on the bubble walls E w . That inconsis- 
tency indicates that in these cases backreaction cannot be neglected. We will comment and 
expand on this issue in section 13.41 

3.2 Fermions 

Turning now to fermionic particle production, in the presence of a tree-level Yukawa coupling 
between the Higgs and the fermions XfH f f, the squared decay amplitude reads 

\M{h^lf)\ 2 = 2\){p 2 -Am)) (3.36) 
which, in the case of SM fermions, leads directly to 



Im 



3 

~ 1 771 ^ f 4 771^ \ ^ 

f{2){x) ]r^k x \ l ~^ L ) e ^- 4m /) (3 - 37) 

The production of (SM) fermions will then be enhanced with respect to the one of Higgs- 
portal S— scalars (specially in the limit of very elastic collisions, see Figure [5]) due to the 
extra factor (x -4m^) in ( 13.371) . Scenarios where the fermionic particle production might 
be important include (apart from the SM itself) the MSSM and its various extensions, due 
to the tree-level coupling between Higgses, Higgsinos and Gauginofl 

In the absence of a direct coupling, the interaction between the Higgs and the fermions 
will occur via an effective operator. This is the case for the so-called fermionic Higgs-portal: 

-AC f = m f Jf + ^\H\ 2 Jf (3.38) 

However, since bubble collisions may excite very massive Higgs field modes (p 2 3> Tew)' 
particle production in this case may be sensitive to the UV completion of the Higgs-portal 
effective theory, making it unreliable to compute the particle production in the fermionic 
Higgs-portal via f)3.38p . Here we consider a simple UV completion for the fermionic Higgs- 
portal, and compute the particle production in this case. We add a singlet scalar field 5* as 
a mediator between the Higgs field and the fermion /, the relevant part of the lagrangian 
being 

- AC f = ^ S 2 + ^ \H\ 2 S 2 + ii s \H\ 2 S + m f Jf + X f Sjf (3.39) 

For simplicity, we will avoid a vev for S (it can be done through the addition of a linear 
term for S in (I3.39P ). For fi s ^ the effective fermionic Higgs-portal operator \H\ 2 ff will 
be generated at tree-level. The squared decay amplitude for h — > f f will then be 

A 2 u 2 v 2 

\M(h->J f)f = 2 ^4-^ (p 2 -4m 2 f ) (3.40) 

with 



4 In particular, the production of neutralino dark matter might have an impact on the subsequent evolution 
of the Universe. 
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Figure 5: Generated comoving energy density T in fermions / (normalized to the observed dark 
matter comoving energy density) as a function of the fermion mass mj in the perfectly elastic 
collision limit (LEFT) and totally inelastic collision limit (RIGHT) for j w = 10 8 , l w = 15/Tew 
and Tew = 100 GeV. Red lines: production in the presence of a direct tree-level Yukawa coupling 
between fermions and Higgs (|3.37p . Blue lines: production for a treedevel effective coupling (|3.42p . 
for fi s = M s = 500 GeV (solid) and 5 TeV (dashed). Yellow lines: production for a ldoop 
effective coupling (|3.43p . The solid-black line corresponds to the observed dark matter comoving 
energy density, and the dashed-black line corresponds to the maximum possible generated comoving 
number density (£ = E w ). 



AM + A* / _ 4mf +^L(l- ) 8(M?-4m}) (3.41) 

16ttM, V Ml K s h> 8tt V M? > K s f> V ; 



leading finally to 



Im 



f (2) (X) 



\2 ,,2 2 
A f fl s U T 



X 



f 4vr ( x -M2) 2 + r2M s 2 



X 



6 ( X - Am)) (3.42) 



When n s = the effective fermionic Higgs-portal operator is not generated at tree-level, 
but rather the decay h —> f f occurs via a finite 1-loop diagram, yielding 



Im 



f (2) (X) 



-- { ^F[ m ),Ml X ] X (l-^ 
(4vr) 5 1 f \ X 



2x1 



6(x-4mJ) (3.43) 



where F [m), M s 2 , x] is a form factor that scales as 
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F [m% Ml X ] — > ^ Log f X » m% M 2 (3.44) 

Fermionic Higgs-portal particle production both in the fi s = and /x s 7^ is shown in 
Figure [3J where it can be clearly seen that the production in the absence of a direct coupling 
between the Higgs and the fermions / differs from what would have been naively obtained 
using (13.380 . As for the case of scalar particle production, under certain circumstances 
the estimate of fermionic particle production neglecting backreaction exceeds the amount of 
energy stored in the bubble walls (S > E w ), and in order to obtain a physically meaningful 
result backreaction should be included (We will expand on this issue in section |3T4"|) . 

3.3 Vector Bosons 

Finally, we study the production of vector boson particles. In the presence of a tree-level cou- 
pling between the Higgs and the vector bosons XyMy hV^V^, the squared decay amplitude 
reads 




(3.45) 



leading to 



{ 3 -ik + Tm)f-*f e( *- 4Mj) (3 - 46) 

Comparing (13.350 . ( 13 . 3 7|) and ( 13.460 we immediately observe the relative efficiency of 
particle production for scalars, fermions and vector bosons. While Im [r^ 2 ^ (%)] scales as x° 
for scalars, and as x f° r fermions, in the case of vector bosons it scales as x 2 i thus greatly 
enhancing production of vector bosons with respect to scalars or fermions for very elastic 
collisions (see Figure [6]) . It is then expected that most of the available energy from the EW 
phase transition will go into and gauge boson production and (possibly) other vector 
bosons coupled at tree-level to the Higgs in extensions of the SmEJ. 

In the absence of a direct coupling, the interaction between the Higgs and the vector 
bosons may occur via an effective operator, as in the so-called vector Higgs-portal [23] : 

- AC V = l -m 2 v VJf* + X v \H\ 2 V^V (3.47) 

However (like for the fermionic Higgs-portal) an analysis of vector boson particle produc- 
tion in the context of the effective theory (13.470 will be unreliable due to very massive Higgs 
field modes (p 2 ^> T| w ) being excited during the bubble collisions. Vector boson particle 
production will then be sensitive to the way in which the effective operator \H\ V^V^ is 
generated. One possible way of generating the effective operator at tree-level, being a 

such as Little Higgs theories or extra-dimensional scenarios with gauge fields living in the bulk. 



Im 



f (2) (X) 
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hidden U(l) gauge field, is by integrating out a £/(l)-charged complex scalar S which has a 

2 

Higgs portal coupling \H\ S*S, the relevant part of the lagrangian then being 

- AC V = -F^F*" - D^S*D^S + V(S) + X hs \H\ 2 S*S (3.48) 

In this scenario, the vector boson acquires a mass via the spontaneous breaking of the 
hidden U (1), through a vev vs for the 5'-scalai@. The squared decay amplitude for h — > 
will then be 
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Figure 6: Generated comoving energy density T in vector bosons (normalized to the observed 
dark matter comoving energy density) as a function of the vector boson mass My in the perfectly 
elastic collision limit (LEFT) and totally inelastic collision limit (RIGHT) for 7^ = 10 s , l w = 
15/Tew and Tew = 100 GeV. Red line: production in the presence of a direct tree- level coupling 
between vector bosons and Higgs (|3.46p . Blue line: production for a tree-level effective coupling 
(I3.50p . for Xh s = 1 and M s = 500 GeV. The solid-black line corresponds to the observed dark 
matter comoving energy density and the dashed-black line corresponds to the maximum possible 
generated comoving number density {£ = E w ). 
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4 (p2 - M s 2 ) 2 + VI Ml 



with T s being the decay width of S. This leads to 



p 



4 My 



(3.49) 



Im 



f (2) (X) 



J v 




9 ( X - AM V ) (3.50) 



6 This implies that there may have been another phase transition in the early Universe associated with 
the spontaneous breaking of the hidden U(l) gauge symmetry, which we need to require to have happened 
long before the EW phase transition since otherwise the EW phase transition would have been effectively 
multi-field and our present analysis of particle production would be totally unrealistic. 
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Vector boson effective Higgs-portal particle production is shown in Figure El resulting in 
a very suppressed particle production with respect to the case in which the vector bosons and 
the Higgs couple directly at tree-level, specially for very elastic collisions. From Figure |6] it 
is also clear that backreaction is most important for direct vector boson particle production 
(for which the production estimate yields £ 3> E w ). 



3.4 Backreaction and Relative Efficiency 

Clearly, for the present analysis of particle production to be physically meaningful it must be 
assumed that the total energy of the produced particles is less than the energy contained in 
the background field configuration h(z,t). Moreover, when the energy of the produced par- 
ticles starts being comparable to the energy of the background field we expect backreaction 
on h(z,t) due to the particle production to be important. Then, in order for the previous 
analysis to be reliable, it is needed 



£ 
A 



< — = -v T — (3.51) 

x Si O L w 



As it has been shown in the previous section, for fermion or vector boson particle produc- 
tion the previous condition f)3.5ip is not satisfied, and in some cases even £ ^> E w is obtained 
(Figure [6] LEFT), signaling the extreme importance of backreaction in those scenarios. 

Since incorporating backreaction into the present analysis of particle production is ex- 
tremely difficult and lies beyond the scope of this paper, we simply note that the relative 
efficiency in particle production for the different species in the present analysis should be 
roughly correct even when backreaction is important. Then, an estimate of the particle 
production in cases where some of the species are very efficiently produced may be obtained 
just by normalizing the production to the total energy in the bubble walls. For very elastic 
bubble collisions, it has been shown in section 13.31 that production of and gauge 
bosons is extremely efficient, which will then leave very little energy left in the bubble walls 
for producing other particle species. The relative efficiencies (defined as ratios of energy in 
produced particles) of the different species for a perfectly elastic collision, normalized to the 
energy contained in the bubble walls (assuming that most of the available energy goes into 
producing and Z^) is shown in Figure [7J A good estimate of the non-thermally generated 
comoving energy density (per particle species a) in this case may then be given by 



x 20 1 (3 £ 
1 ~ — 

v /7f^7 Mpi T EW H A 



u „ 



^ (3.52) 



The fact that this is a reliable estimate of the particle production efficiency for the case 
of very elastic collisions is due to the high-p 2 modes of the bubble wall carrying almost all 
the energy of the bubble wall. The energy carried by the high-p 2 modes will then mostly 
go into vector boson production (their production efficiency at high p 2 is much larger than 
fermionic or scalar ones), result that holds even without incorporating backreaction into the 
analysis. 

On the other hand, for very inelastic collisions the results from the previous section 
show that particle production is only effective for light particles (Mx rrih/2). Therefore, 
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Figure 7: Efficiency of vector boson (solid lines) and fermion (dashed line) particle production 
(scalars are too inefficiently produced to be shown) for a perfectly elastic collision, normalized to 
the most efficiently produced particles (in this case Wn and Zn) and to the energy contained in the 
bubble walls, for j w = 10 8 (LEFT) and j w = 10 15 (RIGHT), l w = 15/T EW and T EW = 100 GeV. 
The solid-black line corresponds to the observed dark matter comoving energy density (normalized 
to the energy contained in the bubble walls). 



production of and will be very suppressed in this case, along with any other heavy 
particle, and most of the available energy will go into production of SM fermions (mainly 
bottom quarks) and (possibly) new light scalars or fermions with sizable couplings to the 
Higgs. 

4 Non-thermal Multi-TeV WIMP Dark Matter 

In this section we focus on the case of relatively heavy dark matter, Mx > TeV, and explore 
the conditions under which the amount of non-thermally produced heavy dark matter can 
end-up accounting for a sizable part of the observed dark matter relic density (dark matter 
may nevertheless still have a thermal component coming from the usual freeze-out process). 
The first condition is clearly that bubble collisions have to be fairly elastic: it has been shown 
in sections [3] that for very inelastic bubble collisions only light (Mx < "Ti/i/2) particles are 
efficiently produced, while heavy particle production is extremely suppressed. Since fast 
thermalization of light species after the EW phase transition seems unavoidable^, for very 
inelastic bubble collisions either dark matter is not efficiently produced or it thermalizes 

1 Dark matter may be coupled to the Higgs weakly enough as to avoid thermalization, however in that 
case we find it is not produced in sufficient quantities to make up for the observed relic abundance. For a 
discussion of the asymmetric dark matter case, see apppendix [Al 



18 



immediately after the end of the EW phase transition, not having any influence on the 
subsequent evolution of the Universe. 

For very elastic bubble collisions, the analysis from sections [3] and 13.41 shows that elec- 
troweak gauge bosons and are most efficiently produced, and the relative production 
efficiency of heavy fermions and scalars is too low (for them to be able to account for a 
sizable part of the observed dark matter relic abundance, see Figure [7]). This leaves heavy 
vector bosons with a direct coupling to the Higgs field as the only viable candidate for 
non-thermally produced dark matter during the EW phase transition. 

In the following we perform an analysis of heavy vector boson dark matter coupled to 
the Higgs, including an overview of thermal freeze-out and direct detection constrains from 
XENON100 [22] (see [231124] for more details), and a comparison between the amount of 
non-thermally produced dark matter and the amount of dark matter produced through 
thermal freeze-out. We also study the evolution of the non-thermally produced dark matter 
component after the EW phase transition. 

4.1 Higgs- Vector Dark Matter Interplay 

Consider a vector boson dark matter candidate with mass My and a tree-level coupling to 
the Higgs [23J[25], 

Ly = X -M 2 V V,V, + XwrhV^V, (4.53) 

This coupling mediates the dark matter annihilation into Standard Model particles, as 
well as the elastic scattering on nucleons relevant for dark matter direct detection. Con- 
cerning the former process, the Higgs boson can mediate annihilation of dark matter into 
electroweak gauge bosons (for heavy dark matter they are the most important annihilation 
channel) through the couplings 

- (2M$rW+W~ + M 2 Z Z,Z,) (4.54) 

The spin-averaged amplitude squared for the annihilation process V^V^ — > W^W~ in the 
limit s ^> ml is given by 

\Mvv^w/z,w,z? - 2 -X v - ^ + 3) (4.55) 

Given (14.551) . the thermally averaged annihilation cross section is given by 

z A 2 f°° fxix — 4) \ 

(™)w^w/z,w/z = 1927rMv V K 2 (zy A dx ^^K^z) [—^- + 3J (4.56) 



where z = My/T, and Ki(z), K 2 (z) are Bessel functions. For z ^> 1, (14.561) reduces to 



A 2 

{°v)yy^ w / Z ,W/Z ~ 167fM 2 ( 4 ' 5 7) 
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The thermal cross section giving rise to the observed value of the relic density {cv) WMAP 
2.6 ■ 10~ 9 GeV~ 2 corresponds, for heavy dark matter My ^> m^ and using (I4.57p . to 
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(4.58) 



WMAP 



Turning now to dark matter direct detection, the spin-averaged amplitude squared for 
Higgs-mediated dark matter elastic scattering on nucleons reads 



\M 



8 Ay/. 



VN^VN\ 



N m Af 



3(t 



m 



2\2 



(My ~ I) 



My 



2m 



N 



16 Xy 



In 



N 



ml 



(4.59) 



Here, m^ m 0.939 GeV is the proton/neutron mass and Jm is the effective Yukawa coupling 
of the Higgs to nucleons which, following [23], we take /jv = 0.326 based on the lattice 



estimate in [26]. In the last step we have taken the limit t <C m 
scattering cross section then reads 
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(4.60) 



On the other hand, the XENON100 bound on the dark matter elastic scattering cross 
section on nucleons, for My > TeV is approximately 



(Jvn^vn < M v ■ 2.2 • 10~ 44 cm 2 (4.61) 

Therefore, (I4.58P and (I4.6ip leave a sizable window in the parameter space (My, Ay) for 
which the dark matter abundance obtained via thermal freeze-out is significantly smaller 
than the observed dark matter relic density, and still the value of Ay (as a function of My) 
is below the XENON100 bound (as shown in Figure [S])- 



4.2 Fate of Non-Thermally Produced Vector Dark Matter 

Given the results from the previous section (summarized in Figure [8]), it is fair to ask if, in 
the region of (My, Ay) parameter space in which the thermal component is not enough to 
account for the observed dark matter relic density, dark matter produced non-thermally at 
the EW phase transition could account for the extra needed amount. Using the results from 
production efficiency of heavy vector boson dark matter obtained in sections 13.31 and 13.4^ we 
show in Figure the value of Ay (as a function of My) for which the amount of non-thermal 
vector boson production equals the observed dark matter relic density (dashed-blue line). 
Then, for values of Ay above the thermal cross section giving the observed relic density 
(af) WMA p, non-thermal production of heavy vector bosons is so efficient as to generate 
amounts of dark matter much larger than the observed dark matter relic density. 

Assuming that at the time of the EW phase transition vector boson dark matter is 
already frozen-out (Tf D ~ My /20 > Tew), we can study the evolution of the non-thermally 
generated dark matter abundance via a simple Boltzmann equation in which the comoving 
dark matter number density Y fulfills Y(z) 3> Yeq(z) (with Yeq(z) being the equilibrium 
comoving number density), yielding 
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Figure 8: The dashed-black line corresponds to the limits on Ay from XENON100 (|4.6ip . The 
solid-black line corresponds to the value of Ay for which the observed DM relic density is obtained 
via thermal freeze-out (14.580 : below it the thermal DM density is larger than the observed DM 
relic density (and thus this region is excluded). Above, the thermal DM density is only a fraction 
of the observed DM relic density, and the red lines show the percentage of relic density accounted 
for by the thermal density. 



with a = (47r 2 i/£ <7*)/45 ~ 2.642 (g* ~ 100 being the number of relativistic degrees of 
freedom in the thermal plasma and £ = 90/(32 7r 3 )), Mpi = 1.2 x 10 19 GeV and y(z) = 
a (a v) Mpi My Y(z). Integration of (14.621) for z > z EW yields 

1111 111 , 

(4.63) 



y{z) v(zew) zew z y(oo) z EW y(z E w) 
Then, given the fact that non-thermal vector boson dark matter production is much 
larger than the observed relic density in the (My, Ay) region of interest, we can take the 
limit y(oo) -C y(z-Ew), obtaining 

y(oo) ~ z EW (4.64) 

From f)4.64p . we immediately obtain that the value of the annihilation cross section that 
will yield the observed dark matter relic density once the non-thermally generated dark 
matter evolves after the EW phase transition is simply given by 

T 

(av) = (av) WMAP — — (4.65) 

-'EW 
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Figure 9: Black lines are the same as in Figure [HJ The dashed-blue line corresponds to the value of 
Ay needed for the non-thermally produced energy density in vector bosons (with a direct coupling 
to the Higgs) to be equal to the DM relic density, for ^ w = 10 8 . The red lines show the values 
of Ay yielding the "non-thermal" cross section (|4.65p (for which the final dark matter abundance, 
taking into account its evolution after non-thermal production, corresponds to the observed dark 
matter relic density) for several values of Tew ■ 



The red lines in Figure [9] show the values of Ay yielding the correct "non-thermal" 
annihilation cross section ( 14.65P for several values of Tew- 

This analysis shows that non-thermal production of multi-TeV vector boson dark matter 
at the EW phase transition (in (My, Ay) parameter space in which the amount of dark 
matter yielded by thermal freeze-out is not enough to account for the observed dark matter 
relic density) is efficient as to generate a dark matter amount much larger than the observed 
relic density. This results in a reactivation of thermalization processes that lead to partial 
wash-out of the non-thermally generated dark matter (wash-out is not complete due to the 
reactivation happening for T < Tew < Tf Q ), meaning that multi-TeV dark matter may have 
a thermal spectrum despite a large fraction of it having been produced non-thermally at 
the EW phase transition. As shown in Figure [91 in the presence of these non-thermally 
produced WIMPs, the relation between mass and coupling giving rise to the observed dark 
matter relic density gets modified with respect to the usual thermal freeze-out scenario, 
leading to better detection prospects in the multi-TeV region for future dark matter direct 
detection experiments. 
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5 Baby-zillas: Super-Heavy Dark Matter from the EW 
Phase Transition 



In this section we study the production of super-heavy dark matter with a mass Mx satisfying 
Mqut ^ M x 3> v ew m the bubble collisions at the end of a very strong EW phase transition. 
We call these dark matter particles baby-zillas because of many similarities (but smaller mass) 
to the WIMP-zilla scenario. 

From Figure El it can be inferred that for 7^ ~ 10 14 — 10 15 non-thermal heavy vector 
boson production in elastic bubble collisions can be so efficient as to generate the observed 
dark matter relic density even for very large dark matter masses My ~ 10 6 — 10 8 GeV and 
perturbative values of the coupling Ay. Using (I3.52p . we plot in Figure ITUl the region in 
parameter space (My, Ay) for which non-thermal production directly yields the observed 
dark matter relic density. 




M y (GeV) 

Figure 10: Region in the (My, Ay) parameter space for which non-thermal production yields the 
observed dark matter relic density for l w = 15/Tew (with Tew = 100 GeV) and j w = 10 14 — 10 15 . 



5.1 Bounds on the Reheating Temperature After Inflation 

A stable particle with mass My ~ 10 5 — 10 8 GeV would yield a much larger relic abundance 
than the observed DM relic density, were it in thermal equilibrium at some stage after 
inflation. For such a massive species, the annihilation cross is always smaller than the one 
needed to yield the observed DM relic density through thermal freeze-out. It is then needed 
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that this particle species never reached thermal equilibrium after the end of inflation. This 
sets an upper bound on the reheating temperature after inflation, specifically Trh < Tf a (with 
Tf being the temperature below which the particle is decoupled from the thermal plasma). 
For a heavy vector boson annihilating into SU(2) gauge bosons (the most important 
annihilation channel in this case) through the Higgs, Tf Q satifies 



where the thermally averaged annihilation cross section (av) is given by (14.561) . In Figure 
ITTI we plot the minimum value of z (corresponding to the maximum allowed value of the 
reheating temperature Trh) as a function of the mass My for the range of Ay values giving 
rise to the observed dark matter relic abundance for 7^ = 10 14 — 10 15 (see Figure fTUI) . We see 
that the upper bound on Trh is relatively insensitive to the precise value of j w , and roughly 
scales as Tgjf ~ My/10. 



Figure 11: Bounds on the Reheating temperature after inflation for the requirement that dark 
matter never reaches thermal equilibrium after inflation, namely Trh < Tf a , as a function of the 
dark matter mass My, and assuming Ay (My) for which non-thermal production yields the observed 
relic abundance (as shown in Figure [TUT) . 



6 Conclusions 

Dark matter may have been efficiently produced at the end of a first order EW phase 
transition if it has a large coupling to the Higgs field. In this paper we investigated the 
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conditions for this non-thermal production mechanism to account for most of dark matter 
in the Universe. We considered scalar, fermion and vector dark matter coupled to the SM 
through the Higgs (either via a direct, tree- level interaction or an effective Higgs-portal 
coupling), and found that production of vector bosons directly coupled to the Higgs is most 
efficient, while for scalar and fermions most of the energy stored in the bubble walls is bound 
to be released into production of SM particles. This analysis singles out vector dark matter 
in the present context. 

For very inelastic bubble collisions only dark with Mx < 100 GeV can be efficiently 
produced, while production of heavier dark matter is extremely suppressed. Unfortunately, 
for a dark matter mass in this range, we did not find a way to avoid subsequent thermalization 
and the wash-out of the non-thermal component, and therefore in this case dark matter 
production at the EW phase transition is irrelevant. The situation is quite different for 
highly elastic bubble collisions. In that case, dark matter with Mx S> 100 GeV can be 
efficiently produced for the so-called runaway bubbles, that expand with a very large 7- 
factor. 

We have identified two scenarios where wash-out of dark matter produced at the EW 
phase transition can be naturally avoided. One has dark matter in the multi-TeV range, 
which makes it possible for non-thermally produced dark matter to remain out of thermal 
equilibrium after the EW phase transition. We determined the region in the parameter 
space of dark matter mass and coupling to the Higgs where the correct relic abundance is 
reproduced. For a given mass, the coupling has to be larger than in the usual thermal freeze- 
out scenario for Higgs portal dark matter, which can be especially relevant for direct detection 
searches, as it opens the possibility of detecting a signal from multi-TeV non-thermal dark 
matter in the near future by XENON100 and LUX experiments. The other scenario is 
baby-zilla dark matter with Mx ~ 10 6 -10 8 GeV. Surprisingly enough, such super-heavy dark 
matter can be produced in important quantities at the end of a strongly first-order EW 
phase transition, provided the dark matter coupling to the Higgs is large, and the 7 factor 
of the bubble walls is near its maximal value of j w ~ 10 15 . In order for the baby-zillas to be 
a viable dark matter candidate, they must have never reached thermal equilibrium, which 
then constrains the reheating temperature after inflation in this scenario. 
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A Asymmetric Dark Matter Production 

We now explore the possibility of asymmetric dark matter production during the EW phase 
transition, together with the viability of this mechanism as a way to avoid wash-out of non- 
thermal production for relatively light dark matter (and any other light species in general). 
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For mult i- component dark matter (X = X a ), an asymmetry in the number densities of 
X a and X a may be generated during the particle production. We will analyze in detail below 
the generation of this asymmetry for scalars (X a = S a ). Then, in section IA.2I we study the 
evolution of the generated asymmetries after the EW phase transition. 

A.l Decay Asymmetries: Producing a Dark Matter Asymmetry 

Let us consider a set of iVj real scalars hi (that includes the field(s) involved in the EW 
phase transition) coupled to a set of N a complex scalars S a via a trilinear interaction. The 
relevant part of the lagrangian is 

- AC = ml S* a S a + C iaP h S* a Sp + V(hi) (A.l) 

where by hermiticity C ia p = C*g a (it follows that C iaa are real, but C ia p with a ^ (3 can 
be complex), and the mass matrix for the scalars S a is taken to be diagonal without loss of 
generality. We also consider a possible term fiijhih^ appearing in V(hi). The lagrangian 
(lA.lj) incorporates a Z 2 symmetry that makes the lightest of the scalars S a stable, which 
may then be a suitable dark matter candidate. In order for an asymmetry in the production 
of S a and S 1 * to be generated, we need a nonzero value for 



\M(hi -+ S* a Sp)\ 2 - \M{hi -> SlS a )\ 2 (A.2) 



At tree level 



j^Tree^fo. _^ c% ^ _ Q*^ =^ \M(hi -)■ 5* Sp)\ 2 — \M(hi ->• Sp S a )\ — \C ia p 



2 



(A.3) 

and there is no asymmetry generated. At 1-loop we include the 1PI diagrams shown in 
Figure [12j Their contribution to the 1-loop decay amplitude is 



It 

3,1,5 

M 1L (h t -> s; s a ) = (c* l5 q a7 q S p i T + ^q a5 q 5p i T ) (a.4) 

J>7>* 

where the integrals It and It correspond to 



(k 2 — m 2 )((k + p) 2 — m 2 )((k — k 2 ) 2 — m 2 ) 



lT = 7T* J d k (k 2 - m 2 )((k + p) 2 - m 2 )((k - k 2 ) 2 - mj) ^'^ 

and can be computed in terms of the usual Passarino-Veltman 3-point scalar loop integral 
Cq. The leading order difference between — > a* (3)\ z and \M(i — > a/3*)\ is due to the 

interference between the tree level and 1-loop decay amplitudes. We obtain 
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\M(hi^ S* a S p )\ 2 - \M(hi^ 313 



Air* 



^ j/iy Im [C* a/3 C jaS Cjsi3\ Im J T + Im [C* aj5 Ci l5 Cj ai C jS p\ Im [J T ] i (A. 6) 



An interplay between "weak" and "strong" phases (complex couplings and imaginary 
part of a 1-loop integral due to particles in the loop going on-shell) is then needed to get 
CP violation in the decay. For Ni — 1 (i = j = 1) both terms on the right-hand side of 
(1A.6j) will vanish for N a < 3 but may be zonzero for N a > 3. For the case of two fields hi 
(Ni = 2), already with two scalars (N a > 2) it is possible to obtain an asymmetry. 




Figure 12: Tree-level and 1PI 1-loop contributions to the decay hi — > S^Sp 



To obtain the total combined production of S a and S 1 * particles, we will just consider the 
tree level contribution to hi — >• S 1 * S a , hi — >■ S 1 * Sp and hi — > Sp S a . We then get 



with 



Im 



f (2) (x) 



|Ciaa| / dU aa + |Cj aj g| / (ill 



(A.7) 
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For the asymmetry in the production of S a and 5** particles we obtain 



(A.8) 
(A.9) 



Im 



f (2) (x) 



Asym 



|M(/ ii ^5 , : < S /3 )| 2 -|M(/ i ,^^^ 



Off 



(A.10) 



with \M{i -> - -> a/5*)| 2 given by fLA6|) . 
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A. 2 Fate of the Generated Asymmetric Abundance 

The asymmetric dark matter production process outlined in the previous section will generi- 
cally result in asymmetries for the comoving number densities for particles and antiparticles 
of the different species A a = Yx a — Yx* 7^ at the end of the EW phase transition. Note 
however that the Z2 symmetry forces the sum of the asymmetries of the different species X a 
to vanish 



5> 







(A.11) 



After the EW phase transition, the comoving number densities for the different species 
Yx a will evolve according to a system of coupled Boltzmann equations. Denoting the sym- 
metric and asymmetric part of the comoving number densities for particles and antiparticles 
of the different species by E a = Y Xa + Y x * 



and A a , we can write 
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(A.12) 
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(A.13) 



where s is the entropy density, z = m L /T {m L is the mass of the lightest species X a ) and 
H(z) is the Hubble parameter. From ( IA.12j) . if the annihilation processes are unsuppressed 
the symmetric comoving number densities for the various species X a will be driven close to 
thermal equilibrium (the small departure from equilibrium being due to the presence of an 
asymmetry A a ) 



(A.14) 



Ideally, in the absence of wash-out of A Q , these processes would delay freeze-out and still 
be active for A a ^> H^ q , annihilating away the symmetric part of the number density and 
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leading to S Q — >• A Q . However, the various processes entering (1A.13j) will tend to erase the 
asymmetries A a . In particular, the process X a + SM — > Xp + SM (responsible for kinetic 
equilibrium among the different species X a ) and the decay process X a — » Xp + SM, will, if 
active, wash-out the asymmetries very rapidly, being also quite insensitive to the temperature 
of the Universe after the EW phase transition (the annihilation process X a + Xp — > SM 
also washes-out the asymmetries, but is suppressed for low Tew)- If these processes are in 
equilibrium after the EW phase transition, they will very rapidly drive 

which, together with flA.llj) . leads to A Q — > 0. While it is possible for the decay process X a — > 
X/3 + SM to be suppressed if the different species are quite degenerate, this automatically 
results in unsuppressed kinetic equilibrium. This seems to rule-out asymmetric production 
as a viable mechanism to avoid complete wash-out of the amount of light dark matter non- 
thermally produced during the EW phase transition, the reason being that it is not possible 
(at least in this simple scenario) to keep all the processes that tend to erase the asymmetries 
A a out of equilibrium after the EW phase transition, while having an efficient particle 
production at the transition itself. 
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